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Abstract 

We explore the role played by the quantum relative phase in a well-known 
model of atom-field interaction, namely, the Dicke model. We introduce an 
appropriate polar decomposition of the atom-field relative amplitudes that 
leads to a truly Hermitian relative-phase operator, whose eigenstates correctly 
describe the phase properties, as we demonstrate by studying the positive 
operator-valued measure derived from it. We find the probability distribution 
for this relative phase and, by resorting to a numerical procedure, we study 
its time evolution. 
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I. INTRODUCTION 



The problem of the interaction of an atomic system with a radiation field is a keystone 
of quantum optics. Needless to say, it is impossible to obtain exact solutions to this problem 
and some approximations must be used; the most common being that the radiation field is 
quasimonochromatic and its frequency coincides almost exactly with one of the transition 
frequencies of the atoms (supposed identical and with no direct interaction between them). 

The two-level atom is the natural consequence of this hypothesis Such an object 
is an important tool because it allows us to describe the matter-field interaction in a very 
simple way, and the results constitute a first step to deal with more realistic situations that 
could include losses, broadening of the atomic lines, etc. To put it bluntly, we can replace 
the whole atomic system by an effective two-level system that accounts for all the relevant 
details of the interaction. The Dicke model ||[|, describing the interaction of A identical 
two-level atoms with a single-mode field in a perfect cavity, is perhaps the archetype of this 
situation. 

In the semiclassical version of this Dicke model, correlations are safely ignored and the 
field is interpreted to be a purely classical electric field [|]-|6|]. Such an approximation has 
proven to be very successful and has the virtue of reducing the problem to the exclusive 
knowledge of the atomic dynamics, which is studied in terms of the inversion and the com- 
ponents of the atomic dipole in-phase and in-quadrature with the field (i.e., the Bloch vec- 
tor). These quadratures are the dispersive and absorptive components of the dipole moment 
effective in coupling to the field. 

For some phenomena, such as spontaneous emission by a fully excited atomic system, 
the quantization of the field is required. Then, one must take care also of the evolution of 
the field amplitudes, but the atomic dynamics is still explained in terms of inversion and 
dipole quadratures. 

Classically, the interaction of matter with light is usually treated within the framework 
of the Lorentz model |7J, which assumes that each electron- ion pair behaves as a single 
harmonic oscillator that couples to the field through its electric dipole moment. In spite 
of its simplicity, it is extraordinarily helpful in developing an intuitive feeling for the phys- 
ical mechanisms involved. Although the dynamics of this model is sometimes expressed in 
quadrature components JIJ, the natural way of understanding its resonant behavior is in 
terms of the relative phase between the field and the atomic dipole ||. 

While the quantum quadratures are well known, and so are the associated eigenstates, the 
operator for this relative phase has resisted a quantum description. At this respect, we think 
that, in spite of its maturity and success, the Dicke model is apparently incomplete since 
it lacks a satisfactory description in terms of this relative phase, indispensable to compare 
with the classical world. The main goal of this paper is precisely the general description of 
that variable. 

When focusing on the relative phase between two subsystems, we think the best way to 
proceed, much in the spirit of our previous work in the subject M, is to try a polar decom- 
position of the quantum amplitudes, which parallels as much as possible the corresponding 
classical factorization. For the relative phase between two harmonic oscillators, this is a 



quite straightforward procedure and leads to a unitary solution |T(J 



For the Dicke model, this polar decomposition seems to be more involved, mainly because, 
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unlike for the case of two harmonic oscillators, the Hamiltonian cannot be cast in terms 
of su(2) operators, but rather in terms of some polynomial deformation of su(2). These 



nonlinear algebras have been examined very recently in quite different physical contexts [ITT 



and, by exploiting these results, it is possible to perform such a decomposition in an elegant 
way, obtaining a bona fide Hermitian operator representing the relative phase we wish to 
examine. 

In this paper we use this operator to introduce the associated probability distribution 
and, then, the most relevant dynamical features of the Dicke model can be easily explained 
using this relevant variable. 



II. QUANTUM DYNAMICS OF THE DICKE MODEL 

The Dicke model describes the interaction of a collection of A identical two-level atoms 
with a quantum single-mode field in a lossless cavity. The spatial dimensions of the atomic 
system are smaller than the wavelength of the field, so all the atoms feel the same field. 
However, the model neglects the dipole-dipole interaction between atoms (i.e., their wave- 
functions do not overlap in the evolution). 

The Hamiltonian for this model in the electric-dipole and rotating-wave approximations 
reads as (in units % = 1) 

H = H + H int , (2.1) 

with 

w f N, 

(2.2) 

A S 3 + g (a f 5_ + aS + ) . 



N = a) a + S 3 (2.3) 

is the excitation number operator, g is the coupling constant (which in this approximation 
is the same for all the atoms and can be chosen as real), and A = cu f — u a is the detuning 
between the atomic and field frequencies. The field mode is described by the annihilation 
and creation operators a and a', while the collective atomic operators are defined by 

S±,3 = £4,3, (2.4) 

and obey the commutation relations 

[S 3 ,S ± } = ±S ± , [S + ,S^] = 2S 3 . (2.5) 

Since all the atoms have the same coupling constant, we need to consider only symmet- 
rical atomic states. Then, let us introduce the atomic Dicke states as 0] 



H 

Hint 



Here 
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|M)a = V ~A\ ^ l+>ii • • • l+U 



3M+1 



JA' 



(2.6) 



where |±)j are the eigenstates of the jth atom and the sum runs over all possible manners 
of choosing M indistinguishable atoms from the group of A atoms. In the space spanned by 
these Dicke states, the action of the collective atomic operators is 



S+ |M) £ 

S- |M) £ 
S3 |M) £ 



(M + 1)(A-M) |M + 1) £ 



M(A-M + \) \M- 1) £ 
(M-A/2) |M) a , 



(2.7) 



where the label M (0 < M < A) denotes the number of excited atoms and —A/2 represents 
the bottom energy level of the atomic system. Therefore, the collective atomic operators 
form a {A + l)-dimensional representation of the algebra su(2) corresponding to a spin A/2. 
The case of a single resonant atom (A = 1) corresponds to the well-known Jaynes-Cummings 
model. 

For simplicity we shall restrict henceforth our attention to case of exact resonance be- 
tween the atomic and the field frequency u a = u>f = uj. Since the field mode is described in 
the usual Fock space |n)f, the natural bare basis for the total system is \n, M) = |n)f <g) |M) a . 
However, it is straightforward to check that 



[Ho, Hi n t] = 0, 



(2.8) 



so both are constants of motion. The Hamiltonian H (or, equivalently, the excitation 
number N) determines the total energy stored by the radiation field and the atomic system, 
which is conserved by the interaction. This means that the appearance of M excited atoms 
requires the annihilation of M photons. This allows us to factor out exp(—iH Q t) from the 
evolution operator and drop it altogether. Hence, we can relabel the total basis as 



\N — M,M) = \N — M) t ® |M) t 



(2.9) 



In such a basis, the interaction Hamiltonian, for a fixed value of N, is represented by the 
tridiagonal matrix 



H 



(N) 
int 



/ h . 

h hx 

hi h 2 

V ; ; ; ••• 



\ 



(2.10) 



with 



h 



M 



yJ{M+l){N -M){A-M). 



(2.11) 



The dimension of this matrix depends on whether A > N or A < N, which are situations 
essentially different and must be handled separately. 

Let us assume that A > N and initially all the atoms are unexcited. Then, M = 
and the conservation of the number of excitations implies that only the states (|2.9| ) with 
< M < N take part in the dynamics. Thus, the dimension of the subspace is N + 1. 
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On the contrary, when A < N the number of initial photons is greater than the number 
of atoms and only the states ( j2.9|) with < M < A are involved in the evolution. The 
dimension is now A + 1 

It is easy to check that, due to the properties of the tridiagonal matrices, the eigenvalues 
are distributed symmetrically with respect to zero, with one eigenvalue equal to zero if there 
are an odd number of them fll2| . 

To find the state evolution we shall need the following matrix elements of the evolution 
operator 

CZ> M (t) = (N-M',M'\exp[-iH^t)\N-M,M), (2.12) 
which can be written as 

CZ> M (t) = E UMjUljeM-ief^}, (2.13) 
j=o 

where U is the unitary matrix that diagonalizes the Hamiltonian and £j are the corre- 
sponding eigenvalues. In what follows we shall use the convention of denoting the dimension 
of the Hamiltonian matrix H-£P by V + 1, that is, 

V = min(N,A). (2.14) 

Now, let us assume that the initial field is taken to be in a coherent state \a){ and that 



the atomic state is initially prepared in an atomic coherent state |0 a []T3l , |T4 



i.e. 



|*(0)) = |a) f ®|C> a , (2.15) 

where 

|«) f = £Q„ |n) f , (2.16) 

n 

Q n being the Poissonian weighting factor of the coherent state (with zero phase) with mean 
number of photons n 



/ n n 

Q n = Je---; (2.17) 



and 



K)a " (1 + \(\r /2 ir ^M!(/-M)! CM|M>a 

A 

= J2 A m |M) a , (2.18) 

M=0 

where the parameter ( is normally rewritten in terms of the spherical angles as 

C = -tan(i?/2)e _< *'. (2.19) 
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In other words, the initial state can be rewritten, taking into account ( |2.9|) , as 

l*(0)> = E Qn-mA m \N-M,M). (2.20) 

N,M 

With this initial condition the resulting state can be recast as 
\V(t)}=exp(-iH int t)\V(0)} 

oo V 

= E E Q N -MA M C M , M (t)\N-M\M'). (2.21) 

N=0M',M=0 

If the initial state is not of the same form, but it has a decomposition with different ampli- 
tudes Am or Q n , equation (|2.21|) is still valid when the appropriate coefficients are taken. 



III. RELATIVE-PHASE OPERATOR FOR THE DICKE MODEL 

In the spirit of our previous work on the relative phase for the Jaynes-Cummings 
model [15|], we shall describe the atom- field relative phase in terms of a polar decompo- 



sition of the complex amplitudes. To this end, let us introduce the operators 

X + = aS + , X- = a'S-, 

(3.1) 

-^3 = S3. 

These operators maintain the first commutation relation of su(2) in ( |2.5|) , [X3, Jf-J = ±X±, 
but the second one is modified in the following way: 

[X.,X+] = P(X 3 ), (3.2) 

where P{X 3 ) represents a second-order polynomial function of the operator X 3 . This is 
a typical example of the so-called polynomial deformations of the algebra su(2). Without 
embarking us in mathematical subtleties, the essential point for our purposes here is that one 
can develop a theory in a very close analogy with the standard su(2) algebra. In particular, 
it is clear that the state \N, 0) plays the role of a vacuum state, since 

X_|iV,0)=0. (3.3) 

Then, we can construct invariant subspaces, as in the usual theory of angular momentum, 
by 

\N-M,M) = X* f |iV,0), (3.4) 

where M is a normalization constant. One can check that 

X^ +1 |iV,0) = 0, (3.5) 
confirming that the number of accessible states is T> + 1. 
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In consequence, the whole space of the system can be split as the direct sum 7i = 
©^ =0 7ijv of subspaces invariant under the action of the operators (X + , X_, X 3 ) and each 
one of them having a fixed number of excitations. These independent subspaces do not 
overlap in the evolution, in such a way that if the initial state belongs to one of them, it will 
remain in that subspace for all the evolution. 

In each one of these invariant subspaces the operator X3 is diagonal, while X + and X_ 
are ladder operators represented by finite-dimensional matrices. This suggests to introduce 
a polar decomposition in the form 

X+ = \J X+X— E 

(3.6) 

X_ = E^X + X_, 

where the radial operator ^X + X_ is diagonal in the basis \N — M, M): 

^Jx + X_\N — M, M) = yfM(N - M + 1)(A - M + 1)\N - M,M), (3.7) 

and 

[X 3 ,E]=E. (3.8) 

We can guarantee now that the operator E = exp(i<£>), representing the exponential of 
the relative phase, is unitary and commutes with the excitation number 

EE* = E^E = 7, 

(3.9) 

[E, N] = 0. 

Thus, we may rather study its restriction to each invariant subspace Hat, we shall denote by 
E^ N \ It is easy to check that the action of the operator E^ in each subspace is given by 

E (N) \N - M,M) = \N- (M+ 1),M + 1), 

(3.10) 

E^\N — M, M) — \N — (M — 1), M - 1). 

Obviously, the action of E^ and E^ N ^ becomes undefined on the marginal states |iV — 
Z>, V) and \N, 0). Therefore, it is necessary to add some conventions for closing the actions 
of these operators on the subspace H-n- By analogy once again with the usual su(2) algebra, 
we shall use standard cyclic conditions and impose (up to global phase factors) 

E {N) \N-V,V) = \N,0), 

(3.11) 

E {N) ^\N,0) = \N-V,V). 
With these conditions, the operator E^ can be expressed as 
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V 



E {N) = j2 \N - (M + 1),M + 1)(N - M,M\ 

M=0 

+ e ^+i)^ ) |iV ) 0)(iV-P,P|, (3.12) 

being an arbitrary phase. Note that the crucial extra term in this equation, which 
establishes the unitarity of E^ N \ is precisely based on the finite number of states. Therefore, 
in each invariant subspace Ti^ there are T> + 1 orthonormal states satisfying 

EW|0W) = e^ ) |0^ ) ), (3-13) 
with r = 0, . . . , V. These states can be expressed as 

\4 N) ) = -7^=j £ e iM ^ N) \N-M,M), (3.14) 
and, by taking the same 2tt window in each subspace, we have 

^-00 + ^, (3.15) 



r 



and 0o is a fiducial or reference phase shift that can be arbitrarily chosen. The expression 
for E on the whole space is 



£=££ w = ££l0fV^ <#°l, (3-16) 



N=0 N=0 r=0 

and, since E is unitary, it defines a Hermitian relative-phase operator 

oo oo T> 

$ = £ $ (A ° = £ £ \4 N) ) 4 N) (4 N) \, (3.i7) 

JV=0 N=0 r=0 

that, obviously, has discrete eigenvalues. In the limit T> ^> 1, this spectrum becomes dense, 
as it might be expected. But, on the opposite limit, one may be surprised to find that the 
state 1 0, 0) is a relative-phase eigenstate (with arbitrary eigenvalue 0o)- While this may 
provide a convincing argument that the theory is unreasonable, we think that is not the 
case. The value of 0o will not lead to any contradictions, because any choice will lead to a 
consistent theory. Our choice of this parameter says nothing about Nature, it only makes a 
statement about our individual preference fT6|j . 

Note as well, that the relative-phase eigenstates are maximally entangled states. This 
has the consequence that the relative-phase operator has no classical correspondence in the 
general case, not even for highly excited states. 



IV. ATOM-FIELD RELATIVE PHASE IN TERMS OF ABSOLUTE PHASES 

In the previous Section we have shown a clear way of obtaining ab initio the atom-field 
relative phase. In spite of this, one could still insist in describing this variable in terms of the 
absolute phases of each subsystem. One must start then from previous descriptions of the 



S 



field and atomic phases and manage them until getting the probability distribution for their 
difference. The goal of this section is to show that this way of proceeding leads naturally 
to a positive operator- valued measure (POVM) []I7|,|n|, and how such a POVM is precisely 



generated by the eigenstates of the relative-phase operator. 

To this end, we shall adopt the elegant axiomatic approach developed by Leonhardt et 



al |L9| to describe the phase properties of both subsystems we are dealing with. To make 
the discussion as selfcontained as possible, we first briefly recall the essential ingredients of 
the general formalism. 

When dealing with generic angle-action variables, one imposes that the complex expo- 
nential of the angle (denoted by E) and the action variable (denoted by L z ) satisfy [compare 



with Eq. 

[L z , E] = E. (4.1) 

If E were unitary, its action on the basis of eigenstates of L z (denoted by \m)) will be as a 
ladder operator 

E\m) = \m + l). (4.2) 
The eigenstates of E (denoted by |0)) provide then an adequate description of the quantum 



angle pD |. However, realistic measurements always involve extra noise beyond that due to 
the intrinsic quantum fluctuations described by canonical conjugation and it is essential to 
extend the canonical formalism by including fuzzy or unsharp generalizations of the ideal 
description provided by E ||21|| . To this end we shall use POVMs, that are a set of linear 
operators A(0) furnishing the correct probabilities in any measurement process through the 
fundamental postulate that 

P{6) = Tr[pA(0)]. (4.3) 
The reality, positiveness, and normalization of P(9) impose 

A f (0) = A(0), A(0) > 0, ( dO A(0) = /, (4.4) 

but, in general, A(0) are not orthogonal projectors like in the standard measurements de- 
scribed by selfadjoint operators. 

In addition to these basic statistical conditions, some other requirements must be imposed 
to ensure that A(0) provides a meaningful description of the angle as a canonically conjugate 
variable with respect L z (even in the sense of a weak Weyl relation ||22|| ). Then, we require 



e W Lz A(0) e~ W Lz = A(0 + 0'), (4.5) 

which reflects nothing but the basic feature that an angle shifter is an angle-distribution 
shifter. This condition restricts the form of the POVM to 

1 oo 

A W = n- E Km e l(m - n)e \m)(n\. (4.6) 

^ n,m=0 
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We must take also into account that a shift in L z should not change the angle distribution. 
A shift in L z is expressed by the operator E, since it shifts the distribution of L z by one 
step. Therefore, we require as well 

E A(0) E ] = A(0), (4.7) 

which, loosely speaking, is the physical translation of the fact that angle should be comple- 
mentary to the action variable. This implies the invariance 

that allows us to recast Eq. ( |4.6| ) as 

1 oo 

A(0) = — E^e-**^, (4.9) 

Z7T — 

v 

while conditions (|4.4|) read now as now 

\K\ < l, K = b u . (4.io) 

Expressing E in terms of its eigenvectors \9), we finally arrive at the general form of a 
POVM describing the angle variable and fulfilling the natural requirements (|4.5|) and Q4.71) : 

A(6) = [ d9' B{9') \6 + 6')(6 + 6'\, (4.11) 

J2n 



where 



1 oo 

B(e) = -J2b u e^. (4.12) 



This convolution shows that this effectively represents a noisy measurement, the function 



B{9) giving the resolution provided by this POVM [20 . 

Let us now focus on the phase properties of our two subsystems. Concerning the field 
phase, the question has attracted the attention of physicists for almost as long as quantum 



mechanics has existed as a physical theory (for recent reviews see Ref. |p3|j ). Nowadays, 



it seems indisputable that an operator representing the phase of a single-mode field in 



a infinite-dimensional Hilbert space cannot exist [24] and the proper way to face up the 
problem involves the use of a relative-state formalism. In spite of this serious drawback, 
a variety of solutions have been proposed to circumvent the difficulties. Virtually all of 
them can be formulated within the POVM formalism discussed before, but with the role of 
L z being played by the number operator a^a. Then a number-shifter is expressed by the 
Susskind-Glogower phase operator |25| (note that we are not concerned about the problems 



of E as a phase operator here, we only use the number-shifter property) and the phase states 



arc 



2tt n=0 



E e f Hf- ( 4 - 13 ) 
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On the other hand, the difficulties with hermiticity that phase operators encounter in 
the case of a single-mode field disappear for a two-level system. In general, for the group 
SU(2) it is possible, by working in the standard (A + l)-dimensional Hilbert space associated 
with the usual angular-momentum operators, to find a truly phase operator from a polar 
decomposition of the amplitudes S± [p6| , p?|1 . The procedure is quite similar to that followed 
in Sec. Ill for the relative phase and, perhaps, the most striking consequence of this approach 
is that the atomic-dipole phase can take only A + 1 different values, due to the dimension 
of the atomic-state space. 

Because of this particular behavior, one may think rather preferable to describe the dipole 
phase by a POVM taking continuous values in a 27r interval, even though this cannot lead to 
a standard operator description. To this end, it suffices to note that the general properties 
( |4.5| ) and (|4.7|) still hold, but the role of L z is played now by S3, and the "Susskind-Glogower" 
eigenstates of the shifter are now 

l^> = 4= E e 4Mea |M) a . (4.14) 
v 2vr M=0 

The joint probability distribution for atomic and field phases can be defined in a very 
natural way as 

P(0 a ,0 f ) =Tr[pA(0 a ,0 f )], (4.15) 

with 

A(0 a ,0 f ) = A a (# a )® A f (£ f ). (4.16) 

Our remaining task is to consistently derive a POVM for the relative phase <fi = 9 & — 6f 
from these expressions. This goal can be achieved in many ways |[28|| : for example, one 



can perform a change of variables to express A(0 a , 9{) in terms of the phase sum and phase 
difference and then remove the phase-sum dependence by simple integration ||20|| . Another 
possibility is to directly define the probability distribution for the relative phase as 

P(4>)= f d9 P(9,9 + <f>) = Tr[pA(0)], (4.17) 

J2tv 



where 



A(0) = / d9 A(9, 9 + </>). (4.18) 

J2n 



For our purposes here, it is sufficient to note that we must get the same values for 
any periodic function of the relative phase whether we use the variable or (9 a ,,9{). In 
consequence, we can impose that 

/ d<f) P(<f>) j"* = I d9 a d9 { P{9 a ,9 { ) e iv ^- di \ (4.19) 

J2-K J2-K 

To proceed with we need to take a definite choice for the corresponding POVMs. Con- 
cerning the field phase, we recall that the Pegg-Barnett formalism and many others embody- 
ing the concept of phase as an observable canonically conjugate to photon number, lead to 
the POVM induced by the Susskind-Glogower phase states, namely fl5| 
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A f (fl f ) = \9 { )(6 { 



(4.20) 



Motivated by this choice, we can use for the atomic phase the finite-dimensional translation 
of this POVM; i.e., 

A a (0 a ) = |0 a )<0 a |, (4-21) 
. A simple calculation shows then that 



A(0)= E I0 (AO ><0 (AO I, (4-22) 

N=0 

where 

1 v 

= — E e iM< ^ | N — M, M). (4.23) 

Therefore, we conclude that the POVM generated by the eigenstates of our relative phase 
operator is just the induced by other absolute-phase approaches (such as Susskind-Glogower 
or Pegg-Barnett), when cast to the appropriate 2tt range. This is, in our view, another 
confirmation that the theory proposed works correctly. 



with \9 a ) given in Eq. (4.14 



V. RELATIVE-PHASE DISTRIBUTION FUNCTION 

For any state, the information one can reap using a measurement of some observable is 
given by the statistical distribution of the measurement outcomes. For the relative phase, 
it seems natural to define the probability distribution function of a state described by the 
density matrix p as 

P(V,0 r ,t) = (0^|p(t)|0W). (5.1) 

However, for physical states |2£j (i.e., states for which finite moments of the number oper- 
ator are bounded) this expression will converge to a simpler form involving a continuous 
probability density we shall write as 

P(V,0,t) = (0W|p(t)|0W), (5.2) 

where the vectors |</>^) defined in (|4.23|) lie in the subspace with total number of 
excitations N. In fact, this expression can be interpreted as a joint probability distribution 
for the relative phase and the total number of excitations. From it, we can derive the 
distribution for the relative phase as 

oo 

P(M= £P(JV,0,t), (5.3) 

N=0 

while 

P(N,t) = [ d(j> P(N,(/>,t) (5.4) 

J2-K 
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can be viewed as the probability distribution of having N excitations in the system. These 
factorizations are an obvious consequence of the fact that the relative phase and the excita- 
tion number are compatible. 

For the general initial state of Eq. (|2.20 ) and the evolution given by Eq. (|2.21 ) we have 



P(iV,0,t) = |(0W|vI/(t))| 2 , 
which, through direct calculation, gives 



P(N,<f>,t) 



1 

2^ 



v 

E Q 

M'M=0 



(5.5) 



(5.6) 



and then we arrive at the total relative-phase probability distribution: 

2 



t) 



2tt 



E 

N=0 



V 



E Qn-mAm C^f, M (t)e 



iM'<t> 



M'M=0 



(5.7) 



This is our basic and compact result which we use to analyze the evolution of the phase 
properties of the Dicke model. 

In Fig. 1 we have numerically evaluated this distribution P{4>, t) as a function of <fi and 
the rescaled adimensional time gt, for the case when all the atoms are initially unexcited 
and the field is in a coherent state with various values of the mean number of photons n. In 
all the cases, when r = we have that C^, M (0) = Sm'm and therefore 



p(<M = o) 



2tt 



E 

7V=0 



V 



E Qn-mAm e 

Af=0 



iM(j> 



(5.8) 



In particular, when all the atoms are initially unexcited only the coefficient Aq survives and 
the previous expression reduces to 



P{<f>,t 



(5.9) 



This flat distribution reflects the fact that the random phase of the dipole in such states 
induces a uniform distribution centred at 0o- At this respect, it is interesting to notice that 
classically the Lorentz model at resonance predicts for the relative phase values of ±7r/2. 
It turns out that this is also a possible choice to fix the reference phase 0o i n the quantum 
description. For simplicity, in all the graphics we have chosen <fto as the origin 0. 

Two quite different behaviors are evident from these graphics. The first occurs in the 
weak-field region f3~0H32"f , when the number of excitations in the system is much smaller 
than the number of the atoms, iV C A. If, for simplicity, we assume that all the atoms 
are unexcited and the average number of photons in the initially coherent field is small, say 
n ~ 1, then we can retain only the dominant terms in Eq. ( |5.7| ), getting 



2vr 



{1 + n[\C^ + \C^ + 2Re(C^f 0]}e" n . 



(5.10) 



We see that, due to the periodic temporal dependence of the terms C^, M (t), this distribution 
is oscillatory for all times, which is corroborated numerically in Fig. l.a. 
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gion 



The second (and perhaps more interesting) case corresponds to the strong-field re- 
-35], when the initial number of photons is much larger than the number of atoms 
A <C iV. Then, following the ideas of Ref. |33] 
can be approximated, up to order A/^/E, by 



one can show that the coefficients C^, M (t) 



— a M'M\ 



-Vt N t), 



where 



N 



2gJN- A/2 + 1/2, 



(5.11) 



(5.12) 



and dj^, M are Wigner d functions which are defined as the matrix elements for finite 
rotations by operators from SU(2) group representations 



(5.13) 



where M, M' = 0, 1, . . . , A. The point now is that essentially only one subspace of dimension 
A + 1 dominates the dynamics. Moreover, a simple calculation using the explicit form of 
these d functions, gives 



oo 

Z/l N=0 



A\ 



^t {{A-M)\M\ 



tan M (n N t/2)e 



iM{4>-n/2) 



cos 2A (n N t/2) } (5.14) 



where we have assumed that all the atoms are initially unexcited. When A ^> 1 and when 
oscillations are well resolved, one can perform an expansion of the square root getting 




J2®N(t) exp 



A' 



A 



- vr/2 + 5 N f 



(5.15) 



where $at(£) is a function of time of complicated structure that accounts for the collapses and 
revivals and that is of little interest for our purposes here, and 5n = a:Tg[ta.n(Q N t/2)}. Now, 
it is clear that, since 6n takes only the values and it, the previous Gaussian distributions 
tend to have two peaks at <ft — i 7r /2, in agreement with the classical expectations. The 
presence of collapses and revivals are evident in Fig. l.b, which confirms previous numerical 
and analytical evidence. The well-known nearly time-independent behavior in the time 
windows between collapse and revival is also clear. As we can see, the distribution tends to 
be randomized in the evolution, although keeping these two peaks at ±7r/2. 

In the intermediate region, when N ~ A, the behavior is more complex, as shown in 
Fig. l.c, and no analytical approximations are available. 

For the particular case of the Jaynes-Cummings model one can diagonalize exactly the 
Hamiltonian in each subspace 7^^, obtaining the well-known dressed states p7| |, that turn 
to be trapping states PS[ ; i.e., the atomic population (5*3 (t)) remains constant in spite of 
the existence of both the radiation field and atomic transitions p£3f| . These states play a 
fundamental role in that model, so it seems interesting to analyze the corresponding problem 
for the case of the Dicke model. In the strong-field limit one can make the replacement 
a — > a = \ffi e l6 and the interaction Dicke Hamiltonian becomes proportional to the operator 



cl 



JO 



S+ + e~ ie S. 



(5.16) 
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where the phase of the classical field has been chosen to coincide with the phase of the initial 
coherent state of the field. The semiclassical atomic states are defined now as eigenstates of 
H c \ taking this phase as zero: 



2S x \P) a = A P \P) £ 



(5.17) 



with A P = A - 2P and P = 0, 1, 



A. 



Following Ref. f33| , we shall call factorized states those states for which the initial field 
is taken to be in a strong coherent state |ot)f and the atomic system is initially prepared in 
a semiclassical atomic state |P) a - For such states, the total wavefunction of the system can 
be approximately written as a product of its field and atomic parts 



!*(*)) -\P(t)) a ®Ht))i 



(5.18) 



with 



\P{t)) a = exp 



\a(t))f = exp 



. A P (S 3 + A/2) 
-i — =gt 



2Jn - A/2 + 1/2 



\P)-< 



(5.19) 



iApJata - A/2 + 1/2 gt 



\a)i 



and one can verify that they are also (approximately) trapping states. For these states, one 
can find after a simple calculation, 



P(<f>,t) 



A + l 



Ea(M|P) a e 



iM<j> 



M 



(5.20) 



The probability distribution is time independent due to the factorization ( |5.18|) . From the 
arguments in Ref. |J3|, one infers that this factorization holds up to times gt ~ \ffi (which 
can be very long times, in this limit) and with an accuracy in the coefficients of the order 

Of A/y/E. 

Moreover, using the properties of the semiclassical atomic states and assuming A> 1, 
one can replace the sum by an integral, obtaining finally 




(5.21) 



i.e., a Gaussian independent of time. In Fig. 2 we have plotted the probability distribution 
obtained from a numerical computation of Eq. ( |5.7| ), showing this quite remarkable behavior, 
except for the presence of very small (almost inappreciable) oscillations superimposed. 

To gain more physical insight in these behaviors, in Fig. 3 we have plotted the evolution of 
the mean value of (sin $) for various values of N, confirming the previous physical discussion. 

To conclude, let us consider the Dicke model in the large-detuning limit; which is usually 
known as the dispersive limit. More specifically, we are in the case when 



A > gVn + IA. 



(5.22) 
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Then, following the procedure developed in Ref. the interaction Hamiltonian in Eq. (|2.2| 
cn be replaced by the effective Hamiltonian 



H eS = AS 3 + \{2a)a + 1)S 3 + \{C - S 2 3 ), 



where 



2 V2 



A 



A' 



(5.23) 



(5.24) 



The obvious advantage of this Hamiltonian is that it is diagonal and allows for a compact 
analytical expression for the coefficients C M , M (t) as 



Cwm^) = 5m>m exp {-it {A(M - A/2) + A [2(N - M) + 1] (M - A/2) 
+ A [C - (M - A/2) 2 ] }) . 



(5.25) 



When the atoms are initially unexcited or excited (or, more generally, when Am = $mk) 
and for any initial state of field, we have 



P(<f>,t) 



1 

27' 



(5.26) 



for all the times. 

For an arbitrary initial state of the atomic system and the field we get 



1 oo 

W) = -£ 

Z ' H N=Q 



M=0 



with 



/# = 2NMX + [A + \{2A + 1)]M 



3\M Z 



(5.27) 



(5.28) 



Since ( p.23| ) is quadratic in the population inversion operator S 3 and is, therefore, analogous 
to the Hamiltonian quadratic in the number operator of a single-mode field propagating 



through a Kerr medium, one could expect jyj that the evolution of coherent atomic states 
in the dispersive limit of the Dicke model leads to the generation of Schrodinger cat states. 
This superposition reaches the most pure form for initial number field states (in particular, 



the vacuum state minimizes the atomic entropy ||42|). 

The situation with the relative-phase distribution is quite different. It is easy to see, for 
example, that if the field is prepared initially in a number state \k) then Q n = 5k n and the 
relative phase distribution is flat. Nevertheless, for initial atomic and field coherent states 
the relative-phase distribution splits for some special times into several humps. These catlike 
states, according to ( |5.28j ), appear at times r = At = 7r/6 (mod27r). To confirm analytically 
this behavior, we expand Eq. (|5.27|) when initially we have strong coherent states for both 
field and atoms, with n ^> A ^> 1. By replacing once again the sum by integrals, one easily 
gets 



P(0,t = 7T/6A) 




| e _[0_ 7r9iri /(3A)] 2 A/2 + e -[</,+7r-7r 9 i fi /(3A)] 2 yl/2| ^ 



(5.29) 
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where 



(p n = 2n\ + A + X(2A + 1), (5.30) 

and all the phases must be understood mod(27r). The two separated Gaussians indicates 
the presence of two humps and, therefore, the presence of catlike states. To further confirm 
this, in Fig. 4 we have computed numerically the distribution function P(4>,t = 7r/6A) at 
the times predicted by the theory. The graphic clearly demonstrates the presence of the 
two-component state, according to our previous considerations. 

VI. CONCLUSIONS 

In this paper we have investigated an appropriate operator for the quantum description 
of the relative phase in the Dicke model. We have used a proper polar decomposition of the 
corresponding field amplitudes, much in the spirit of our previous work on the subject. This 
polar decomposition has been justified on physical grounds, as well as using the theory of 
polynomial deformations of su(2). 

The eigenvalue spectrum of this operator is discrete, as it happens for the polar decom- 
position corresponding to two field modes. From these eigenstates we have obtained the 
probability distribution for the relative phase and we have studied its time evolution. For 
the weak-field region, the behavior is essentially oscillatory, while for the strong-field region, 
the relative phase tends to be randomized in the evolution, although showing collapses and 
revivals. For both limiting regions, we have developed analytical approximations that have 
allowed us an easy physical interpretation of some remarkable phenomena. 
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FIGURES 

FIG. 1. Gray-level contour plot of the probability distribution P(<f),i) as a function of 4> and 
the rescaled time gt for the case of A = 5 atoms initially unexcited and the field in a coherent state 
with the following values of the mean number of photons: a) n = 1 (weak field), b) n = 50 (strong 
field), and n = 5 (intermediate field). 

FIG. 2. Probability distribution function P(4>, t) as a function of (ft and the rescaled time gt for 
the case of a factorized state with A = 3. The atomic coherent state has -d = tt/2 and ip = and 
the field state has n = 20. 

FIG. 3. Plots of (sin</>) versus for the same values of n as in Fig. 1 (from top to bottom). 

FIG. 4. Probability distribution function P(4>,t) as a function of cf> for the time Xt = tt/6 for 
the case of an atomic coherent state (A = 5) with # = tt/2 and ip = and a field coherent state 
with n = 10. The presence of the two humps corroborate the presence of a catlike state. 
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